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Abstract 

Sympatric cryptic species, characterized by low morphological differentiation, pose 
a challenge to understanding the role of interspecific competition in structuring 
ecological communities. We used traditional (morphological) and novel molecular 
methods of diet analysis to study the diet of two cryptic bat species that are sympatric 
in southern England (Plecotus austriacus and P. auritus) (Fig. 1). Using Roche FLX 
454 (Roche, Basel, CH) high-throughput sequencing (HTS) and uniquely tagged 
generic arthropod primers, we identified 142 prey Molecular Operational Taxo- 
nomic Units (MOTUs) in the diet of the cryptic bats, 60% of which were assigned 
to a likely species or genus. The findings from the molecular study supported the 
results of microscopic analyses in showing that the diets of both species were dom- 
inated by lepidopterans. However, HTS provided a sufficiently high resolution of 
prey identification to determine fine-scale differences in resource use. Although 
both bat species appeared to have a generalist diet, eared-moths from the family 
Noctuidae were the main prey consumed. Interspecific niche overlap was greater 
than expected by chance (Ojk = 0.72, P < 0.001) due to overlap in the consumption 
of the more common prey species. Yet, habitat associations of nongeneralist prey 
species found in the diets corresponded to those of their respective bat predator 
(grasslands for P. austriacus, and woodland for P. auritus). Overlap in common 
dietary resource use combined with differential specialist prey habitat associations 
suggests that habitat partitioning is the primary mechanism of coexistence. The 
performance of HTS is discussed in relation to previous methods of molecular 
and morphological diet analysis. By enabling species-level identification of dietary 
components, the application of DNA sequencing to diet analysis allows a more 
comprehensive comparison of the diet of sympatric cryptic species, and therefore 
can be an important tool for determining fine-scale mechanisms of coexistence. 



Introduction 

Interspecific competition is an important mechanism struc- 
turing ecological communities (Schoener 1983), outweigh- 
ing the effect of positive mutualistic interactions (Alexan- 
drou et al. 2011), and promoting increased speciation rates 
(Thierry et al. 2011). Stable coexistence in the face of inter- 
specific competition is the result of ecological niche differen- 
tiation, whereby phenotypic character divergence promotes 
differential resource use, thus reducing the effect of compe- 



tition for limiting resources (Chesson 2000). There is ample 
evidence to support both ecological character displacement 
(species occupying different niches in sympatry; reviewed 
in Dayan and Simberloff 2005) and limits to the similarity 
of coexisting competitors (e.g. Kingston et al. 2000). More- 
over resource partitioning has been identified as a mecha- 
nism that facilitates coexistence in a variety of animal com- 
munities (reviewed in Schoener 1974). Nevertheless, some 
studies have shown that ecologically similar species can co- 
exist through neutral processes that promote phenotypic 
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character convergence rather than divergence (Leibold and 
McPeek 2006). 

Traditional views of niche-based species coexistence are 
challenged by the presence of sympatric, morphologically 
similar, but genetically isolated, species (cryptic species) that 
do not appear to differ sufficiently in their morphology to 
allow niche differentiation (Wellborn and Cothran 2007). 
However, morphological similarity may not be a sufficient 
indication of ecological similarity, and even small differ- 
ences can influence access to resources (Saunders and Bar- 
clay 1992). Cryptic species appear to be common in bats 
(order Chiroptera) (Clare 2011; Clare et al. 2011a), and in 
particular among insectivorous bats, because differences in 
their sensory abilities, in the form of ultrasonic echoloca- 
tion calls, are not readily distinguished by humans. More- 
over, echolocation gives bats an added dimension for niche 
separation that may be less conserved in its evolution than 
morphology per se (Jones 1997), and therefore may pro- 
mote ecological divergence potentially leading to speciation 
(Kingston and Rossiter 2004). Because foraging habitat use 
and prey selection in insectivorous bats are closely linked to 
wing morphology (Aldridge and Rauthenback 1987; Norberg 
and Rayner 1987) and echolocation call structure (Jones and 
Rydell 2003), even moderate differences in these features may 
allow species to partition resources. 

Despite their morphological similarity, sympatric cryp- 
tic bat species may show pronounced spatial segregation of 
primary foraging habitats. For example, partitioning of for- 
aging habitat, in the absence of apparent differences in wing 
morphology (Jones 1997), was identified among the recently 
separated cryptic Pipistrellus species, whereby P. pygmaeus 
selects riparian habitats, while P. pipistrellus prefers a wide 
range of habitats including deciduous woodlands and pas- 
ture (Davidson-Watts et al. 2006; Nicholls and Racey 2006). 
Differences in echolocation call structure can facilitate niche 
differentiation in sympatric bat species by affecting the detec- 
tion distances of prey of various sizes (Kingston and Rossiter 
2004), as well as contributing to foraging performance un- 
der different levels of habitat clutter (Siemers and Schnitzler 
2004). For bats foraging in the same habitat, sensory adap- 
tations may allow access to different prey, thus reducing re- 
source competition (Siemers and Swift 2006). However, the 
role of interspecific competition in driving trophic resource 
partitioning in bats is still debated (Husar 1976; Arlettaz et al. 
1997; Schoeman and Jacobs 2011). 

The advent of finer resolution methods of diet analysis may 
offer a clearer picture of trophic resource partitioning among 
coexisting species. Indeed some studies that failed to iden- 
tify resource partitioning among morphologically similar bat 
species suggest that finer scale differences in microhabitat 
use or finer resolution diet analysis may reduce the extent 
of resource overlap (Jacobs and Barclay 2009; Schoeman and 
Jacobs 201 1). Contrary to traditional methods of diet analy- 
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sis (e.g., fecal or stomach content analysis), polymerase chain 
reaction (PCR) amplification and DNA sequencing of prey 
remains from predator fecal samples allows identification to 
the species level of both soft- and hard-bodied consumed 
prey (Symondson 2002). Despite the high potential of DNA 
sequencing to aid identification at finer resolution than the 
ordinal level for important insect prey, such as Lepidoptera 
(Whitaker et al. 2009), this technique has only recently been 
applied to study the diet of insectivorous bats. However, 
molecular diet studies to date have focused on individual 
species (Clare et al. 2009, 2011b; Zeale et al. 2011) or on 
sympatric species that are not close relatives (Bohmann et al. 
201 1), rather than on cryptic species, where finer resolution 
of prey types may be fundamental to decipher mechanisms 
of resource partitioning. 

We applied traditional (morphological) and novel next 
generation sequencing approaches to study the diet of two 
sympatric cryptic sister species in southern England, the gray 
and brown long-eared bats, Plecotus austriacus and P. auri- 
tus. Plecotus austriacus (Fig. 1) is one of the rarest mammals 
in Britain, with a prebreeding population estimated at 1000 
individuals, and is restricted to southern England, while P. au- 
ritus is more common and widespread throughout the British 
Isles (Harris et al. 1995). The long-eared bat genus contains 
at least 19 distinct cryptic species, nearly all of which were 
only identified in the past decade based on molecular studies 
(Spitzenberger et al. 2006; Mayer et al. 2007). Though rec- 
ognized as separate species in the 1960s (Corbet 1964), P. 
austriacus and P. auritus show high overlap in most morpho- 
logical characteristics (Ashrafi et al. 2010), including their 
wing morphology (Sevcik 2003) and echolocation call pa- 
rameters (Russo and Jones 2002). Significant differences in 
the shape of the baculum (Dietz et al. 2009) may partly ac- 
count for reproductive isolation between these cryptic species 
(Patterson and Thaeler 1982). The two bats are sympatric in 
southern England and parts of Europe. However, while P. 
austriacus is primarily a southern European species, P. auri- 
tus is abundant in central and northern Europe, but confined 
to mountainous areas with cooler climates in southern Eu- 
rope (Spitzenberger et al. 2006). Similarities in morphology 
and echolocation calls suggest that these species share several 
niche dimensions, and therefore studying their patterns of 
resource use may reveal important fine-scale mechanisms of 
coexistence (Tokeshi 1999). 

Our objectives were to assess the ability of PCR-coupled 
next generation (Roche FLX, Roche, Basel, CH) sequencing 
to identify fine-scale differences in trophic ecology from fe- 
cal samples, and compare the performance of this technique 
with previous molecular (PCR-coupled cloning and fragment 
targeting) and nongenetic methods of diet analysis. We hy- 
pothesized that despite showing few phenotypic differences, 
cryptic species are able to coexist due to differences in their 
ecology. Therefore, we predicted that (1) coexistence among 
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Figure 1. (a) The grey long-eared bat, Plecotus austriacus, and (b) the 
brown long-eared bat, Plecotus auritus. 

P. austriacus and P. auritus is mediated via fine-scale mecha- 
nisms of resource partitioning; (2) the extent of dietary niche 
overlap and diet selectivity are affected by variation in prey 
availability; and ( 3 ) the diet of the two bat species corresponds 
to patterns of foraging habitat selection. 

Methods 
Sample collection 

We collected fecal samples directly from 30 P. austriacus (rep- 
resenting approximately 3% of the prebreeding British pop- 
ulation [Harris et al. 1995]) and 30 P. auritus bats caught at 
sites in southern England (Devon, Isle of Wight, and Somer- 
set). The P. austriacus samples were collected between April 
2009 and October 2010, and the P. auritus samples between 
June and October 2010 (Supporting Information 1). Bats 
were caught under license either in roosts, bat boxes, or in 
woodlands, and held in separate holding bags for a maximum 



O. Razgour etal. 

of 30 min. Feces were immediately stored in 100% ethanol 
in individual tubes, and frozen within 12 h of collection. For 
the morphological diet analysis, feces were collected once a 
month between April and October 2009 from seed collection 
trays placed on the floor of a P. austriacus maternity colony 
in Devon (50°3'N, 3°3'W). 

DNA extraction and PCR amplification 

DNA was extracted from feces using the QIAamp DNA stool 
extraction kit (Qiagen, Valencia, CA) following the modi- 
fications described in Zeale et al. (2011). A short (157 bp 
excluding primer) region of the cytochrome c oxidase sub- 
unit I (COI) mitochondrial DNA gene was subsequently PCR 
amplified from each DNA extract. We used generic arthropod 
primers with a wide taxonomic coverage that includes 13 in- 
sect and arachnid orders commonly found in insectivorous 
bat diets (ZBJ-ArtFlc and ZBJ-ArtR2c; Zeale et al. 2011). 
The primers were modified into 5'-tagged "fusion primers," 
in order to enable Roche FLX sequencing of pools of the am- 
plicons, and subsequent bioinformatic sorting into original 
PCR (Binladen et al. 2007). Post-PCR amplification, the am- 
plicons were purified, quantified by real-time PCR (qPCR) 
(Meyer et al. 2007), and then pooled by species at equimo- 
lar ratio. Each pool of amplicons was deep sequenced on 
one-eighth of a Roche GS-FLX platform using Titanium se- 
quencing chemistry. PCR reaction conditions, cycle program, 
and sequencing procedures followed Bohmann et al. (2011) 
and are provided in Supporting Information 2. 

Analysis of FLX sequencing data 

We used a conservative approach discarding sequences that 
did not match primers exactly, fragments that were too short, 
and all singleton sequences (unique sequences that were only 
present once) to minimize the effect of sequencing errors. Se- 
quences were aligned using the ClustalW algorithm in BioEdit 
(T. Hall, http://www.mbio.ncsu.edu/bioedit/bioedit.html), 
and the nucleotide and amino acid alignment was checked 
by hand with reference to known arthropod sequences. All 
sequences were collapsed into Molecular Operational Tax- 
onomic Units (MOTUs), approximating species, using the 
software jMOTU (Jones et al. 201 1, https://www.nematodes. 
org/bioinformatics/jMOTU/). This approach is useful for 
describing niche breadth and overlap when identification 
of sequences is limited by an incomplete reference dataset 
(Clare et al. 2011b). The MOTU assignment was tested at 
a series of thresholds corresponding to 0-10% sequence di- 
vergence and the resulting assignments were graphed. For 
further analysis, we selected a threshold (2%) based on the 
infliction point of the graph of distribution of MOTUs re- 
covered, that is, the percent of sequence differences where the 
graph reached an asymptote and new MOTUs were not be- 
ing recognized but existing MOTUs had not been collapsed. 
This value is biologically meaningful and is well within ex- 
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pected variation at the species level in surveyed insects (see 
Clare et al. 2011b for a discussion of MOTU divergences). 
Representative sequences from each MOTU (the three 
most common sequences when identical sequences are col- 
lapsed into haplotypes) were compared to reference se- 
quences using the Barcode of Life Data System (BOLD, 
www.barcodinglife.org) identification engine. We used a con- 
servative threshold of >98.5% sequence similarity for species 
identification based on mean sequence divergence values es- 
timated by Zeale et al. (2011) for the amplified COI region. 

Phylogenetic approaches employed by Clare et al. (2009, 
2011b) to compliment distance-based species identification 
were not possible due to the short length of the sequences 
(157 bp). Moreover, short sequence length meant that some 
MOTUs gave >98.5% match to more than one species, some 
of which belonged to different genera or even families. To 
overcome these obstacles, we used the following criteria to 
create three identification confidence levels: 

1. Solid match (>98.5%) to one species — species-level as- 
signment, or match (>98.5%) to more than one species, all 
belonging to the same genus — genus-level assignment. 

2. Match (>98.5%) to more than one species belonging to 
different genera, only one of which was a U.K. species — 
species-level assignment to U.K. species. 

3. Match (> 98%) to several species of different genera within 
the same family or to reference sequences only identified to 
the family-level — family level assignment. 

Diet analysis 

The diet of P. austriacus and P. auritus was quantified based 
on the percent frequency of occurrence (%FO) of prey orders 
(the number of feces containing an order divided by the total 
occurrences of all orders). All other analyses were carried out 
at the MOTU (species) level. 

The extent of dietary specialization and diversity was de- 
termined at the MOTU level using the standardized Levins' 
measure of niche breadth (Equation 1) and Shannon's diver- 
sity index (Equation 2). 

in 



B : 



pi 



Standardized as : B A = 



B - 1 
n- 1 



where B is Levins's measure, p; is the proportion of fecal 
samples in which MOTU i was found, and n is the number 
of possible MOTUs in the diet. 



H' = -J2piLn(Pi), 



(2) 



We used Pianka's (1973) measure of niche overlap (Equa- 
tion 3) to quantify dietary resource overlap at the MOTU 
level between the two bat species and within each bat species 
between the two main collection sites (Devon and Isle of 
Wight). Null models were used to test whether the extent 
of niche overlap is greater than expected by chance, and 
determine the effect of season and sex on dietary resource 
use. We generate 10,000 simulated matrices of randomized 
MOTU diet composition, using the software EcoSim (version 
7; http://grayentsminger.com/ecosim.htm) with Randomisa- 
tion Algorithm 3 , and compared observed and randomly sim- 
ulated extents of niche overlap. Bonferroni corrections were 
applied to retain significance value at P < 0.05, resulting in 
significance level set at P < 0.017. 



O jk = 



(3). 



P 2 ij P 2 ik 



where p; is the proportion of fecal samples in which MOTU 
i was found. 



where P y is the proportion that resource i is of the total 
resources used by species Pa is the proportion that resource 
i is of the total resources used by species k; and n is the total 
number of resource states (total number of MOTUs). 

Dietary prey composition was compared within each 
bat species between seasons using nonparametric tests. We 
compared the proportion of Lepidoptera MOTUs within 
each fecal sample between the spring (April-May), summer 
(June- August), and autumn (Sep-Oct) for P. austriacus, and 
between the summer and autumn for P. auritus. Statistical 
analyses were carried out in PASW Statistics 18. 

Arthropod prey habitat associations and presence in south- 
ern England were determined based on previously published 
information (Heath and Emmet 1983; Waring et al. 2003; 
Chinery 2005; Manley 2008; Chandler 2010). Arthropod 
species commonly found in a range of habitat types were 
classified as habitat generalists, while species that were exclu- 
sively associated with either grassland and open habitats or 
woodland habitats were classified as grassland and woodland 
specialists, respectively. 

Morphological diet analysis 

Primer generality 

In order to test for primer amplification biases against certain 
prey orders, we used traditional microscopic fecal analysis 
methods (Whitaker et al. 2009) to determine prey compo- 
sition in all fecal samples analyzed genetically. We looked 
for the presence of prey orders that were absent from the 
molecular study and for differences in prey order composi- 
tion in terms of percent frequency of occurrence (%FO) in 
diet. 
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Traditional morphological diet analysis 

We compared the results of the high-throughput sequencing 
(HTS) analysis to a previous microscopic fecal analysis study 
of the diet of P. auritus from Somerset, southern England 
(Hollyfield 1993). However, because there were no previous 
studies of the diet of P. austriacus in Britain, we analyzed 
feces collected from a maternity colony in Devon, using tra- 
ditional microscopic fecal analysis methods (Shiel et al. 1997; 
Whitaker et al. 2009). We randomly selected 30 droppings 
for analysis from each monthly collection, and teased them 
apart under a dissection microscope, picking out all recog- 
nizable arthropod fragments for identification with reference 
slides and identification guides (Shiel et al. 1997; M. James 
pers. comm.). The diets of P. austriacus and P. auritus were 
quantified based on percent frequency of occurrence of prey 
orders (%FO). 



Results 

Of the 30 fecal samples collected for each bat species, 28 
of the P. austriacus (93%) and 24 of the P. auritus extrac- 
tions (80%) produced PCR amplicons, and were therefore 
included in subsequent analyses. The sequencing run gen- 
erated 100,508 and 73,606 sequences containing correctly 
barcoded primers for the P. austriacus and P. auritus samples, 
respectively. Subsequent filtering of the data to return only 
high-quality, identifiable arthropod sequences and removal 
of singleton sequences reduced these to 91,476 and 63,715, 
respectively. 

The high-quality arthropod sequences were collapsed into 
142 MOTUs (using a 2% sequence divergence threshold), 
belonging to three species of arachnids, two species of Crus- 
tacea, and at least six orders of insects. We identified 60% 
of the MOTUs to the likely species or genus (69% were 
solid matches — confidence level 1, while 31% matched more 
than one species, only one of which was a U.K. species — 
confidence level 2), and 16% to the family level (Table 1). 
The remaining MOTUs (24%) could not be identified confi- 
dently based on our classification scheme and were therefore 
left as unknown. In Lepidoptera, where the reference dataset 
is relatively complete, 59 discrete MOTUs were identified as 
52 separate species; suggesting that using current taxonomic 
classification in BOLD, MOTUs slightly overestimated species 
richness (by 12%). 

The majority of MOTUs (80%) were consumed by a single 
individual, though the most common prey species in the 
diet of both species, Noctua pronuba, was consumed by 19 
P. austriacus and 13 P. auritus bats. Number of prey taxa 
(MOTUs) per dropping ranged between 1 and 17 (means ± 
SD: P. austriacus, 5.3 ± 4; P. auritus, 4.2 ± 3). 
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Diet composition and specialization 

The diet of P. austriacus contained 90 MOTUs, while that of 
P. auritus contained 68. Analysis carried out at the MOTU 
level showed that both bat species had relatively narrow 
niches (Levins' measure: B A = 0.21, for both species) but 
high dietary diversity (Shannon diversity index: P. austriacus 
H = 4.1, P. auritus H = 3.9). We identified six prey orders in 
the diet of P. austriacus (Lepidoptera, Diptera, Neuroptera, 
Hemiptera, Trichoptera, and Isopoda) and five orders in the 
diet of P. auritus (Lepidoptera, Diptera, Coleoptera, Araneae, 
and Isopoda). The diet of both bat species comprised pri- 
marily of prey from the orders Lepidoptera and Diptera, with 
Lepidoptera accounting for 66.7% of the diet of P. austriacus 
and 64.7% of the diet of P. auritus (Fig. 2). The same prey 
orders were identified at similar proportions in the morpho- 
logical analysis of prey remains in genetically analyzed fecal 
samples, with the exception of Coleoptera, which was only 
found in the morphological analysis of the P. austriacus sam- 
ples. However, it is important to note that Coleoptera was 
only identified in one P. austriacus fecal sample, and com- 
prised less than 5% of the volume of that sample. 

The majority of the lepidopteran species identified in the 
diet of both bat species were eared-moths of the family 
Noctuidae (proportion of Lepidoptera MOTUs: P. ausitra- 
cus 71%, P. auritus 83.3%; Fig. 3). Species of the family Tip- 
ulidae were the most common dipterans identified in the 
diet of P. austriacus (76%), while species of the dipteran 
family Muscidae were more common in the diet of P. au- 
ritus (65%). The most common prey species in the diet of 
P. austriacus were N. pronuba, Autographa gamma, Hepialus 
sp., Apamea monoglypha, and Tipula oleracea, while the most 
common prey species in the diet of P. auritus were N. pronuba, 
Xestia xanthographa, Hepialus sylvina, and Polietes lardarius 
(Table 1). 

The diet composition at the ordinal level of both bat species 
differed between the seasons. Plecotus austriacus consumed 
higher proportion of Lepidoptera in summer than in spring 
(Kruskal-Wallis test: H = 7.72, N = 14, 8, P = 0.002; Fig. 4), 
while P. auritus consumed more Lepidoptera in summer 
than in autumn (Mann-Whitney test: U = 33, N = 10, 
14, P = 0.014; Supporting Information 3). 

Dietary niche overlap 

We found little evidence for dietary resource partitioning 
among these cryptic bat species. Interspecific niche overlap, 
measured based on MOTUs, was significantly higher than 
expected by chance (Ojk = 0.72, P < 0.001) and higher 
than intraspecific niche overlap (P. austriacus: Ojk = 0.63, 
P = 0.009; P. auritus: O jk = 0.31, not significant) (Fig. 5). 
Niche overlap between bat species was greater in the summer 
(O jk = 0.62, P = 0.012) than the autumn (O jk = 0.48, not 
significant). 
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Table 1. List of arthropod prey identified in the feces of 28 Plecotus austriacus and 24 P. auritus from southern England, includ 


ng the number of 


individuals from each bat species that consumed the prey taxa. Confidence levels are based on the BOLD identification system, whereby confidence 


level 1 = solid match to one species or genus (> 


98.5%); level 2 = match to more than one species ( 


>98.5%), only one of which was a U.K. species; 


and level 3 = match > 98% to several species of different genera, or to reference sequences only identified to the family level. 




Order Family 


Species 


P. austriacus 


P. auritus 


Confidence level 


Lepidoptera Aractiidae 


Unknown 


2 


0 


3 




Diaphora mendica 


1 


0 


1 




Seirarctia sp. 


1 


0 


l 




Spilosoma luteum 


3 


0 


i 


Crambidae 


Agriphila tristella 


2 


0 


1 




Chrysoteuchia culmella 


2 


0 


2 




Crambus perlella 


1 


0 


-y 
Z 


Elachistidae 


Depressaria daucella 


1 


0 


1 
I 




Semioscopis sp. 


2 


0 


I 


Gelechiidae 


Carpatolechia decorella 


1 


2 


I 


Geometridae 


Gymnoscelis rufifasciata 


1 


0 


2 


Hepialidae 


Hepialus sp. 


6 


2 


-y 
Z 




Hepialus sylvina 


6 


3 


-1 
1 


Noctuidae 


Unknown 


6 


3 


3 




Acronicta alni 


0 


1 


1 




Agrochola litura 


0 


1 


1 




Agrochola lota 


0 


1 


I 




Agrochola lychnidis 


0 


1 


Z 




Agrochola macilenta 


0 


1 


1 




Agrotis exclamationis 


4 


2 


1 




Agrotis ipsilon 


2 


0 


-y 
Z 




Agrotis segetum 


1 


0 


1 




Agrotis puta 


1 


0 


1 




Allophyes oxyacanthae 


1 


2 


Z 




Apamea crenata 


0 


1 


1 




Apamea epomidion 


0 


1 


I 




Apamea monoglypha 


6 


0 


-y 
Z 




Apamea sp. 


1 


0 


I 




Autographa gamma 


11 


2 


1 




Charanyca trigrammica 


1 


0 


Z 




Conistra sp. 


0 


1 


2 




Diarsia sp. 


0 


2 


z 




Hoplodrina ambigua 


1 


0 


1 




Hydraecia micacea 


1 


0 


z 




Hydraecia sp. 


1 


0 


-1 
1 




Hypena proboscidalis 


0 


1 


I 




Lithophane hepatica 


1 


0 


2 




Mythimna albipuncta 


2 


0 


i 




Mythimna pallens 


0 


1 


z 




Noctua comes 


2 


0 


2 




Noctua pronuba 


19 


13 


1 




Noctua sp. 


1 


1 


z 




wLiuupieuia sp. 


u 


1 

1 


z 




Oligia sp. 


0 


1 


1 




Omphaloscelis lunosa 


1 


1 


1 




Orthosia gothica 


1 


0 


2 




Phlogophora meticulosa 


2 


0 


1 




Scoliopteryx libatrix 


1 


0 


1 




Tholera cespitis 


0 


1 


1 




Tholera decimalis 


2 


0 


1 




Xestia xanthographa 


2 


6 


2 


Oecophoridae 


Diurnea fagella 


0 


1 


2 




Endrosis sarcitrella 


0 


1 


2 



(continued) 
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Table 1. Continued 



Order 


Family 


Species 


P. austriacus 


P. auritus 


Confidence level 


Diptera 


Anisopodidae 


Unknown 


0 


1 


3 




Calliphoridae 


Unknown 


0 


1 


3 






Pollenia sp. 


0 


1 


2 




Muscidae 


Unknown 


0 


2 


3 






Eudasyphora cyanella 


1 


2 


1 






Eudasyphora sp. 


0 


1 


1 






Helina impuncta 


1 


0 


1 






Morellia simplex 


0 


1 


1 






Phaonia subventa 


1 


3 


1 






Polietes lardarius 


0 


3 


1 




Scathophagidae 


Scathophaga stercoraria 


1 


1 


1 




Syrphidae 


Platycheirus sp. 


0 


2 






Tachinidae 


Leskia sp. 


1 


4 


1 




Tipulidae 


Unknown 


1 


0 








Tipula oleracea 


8 


1 


1 






Tipula sp. 


1 


0 


1 


Neuroptera 


Chrysopidae 


Unknown 


2 


0 




Coleoptera 


Dermestidae 


Anthrenus fuscus 


0 


1 


1 


Worn i n+cira 
nfcM 1 lip LfcM d 


IVI 1 1 lUdc 


U 1 1 l\l IUWI 1 




u 








Phytocoris tiliae 


1 


0 




Trichoptera 


Limnephilidae 


Stenophylax sp. 


1 


0 




Isopoda 


Porcellionidae 


Porcellio scaber 


1 


1 








Porcellio sp. 


0 


1 




Araneae 


Anyphaenidae 


Unknown 


0 


1 








Anyphaena accentuata 


0 


1 






Tetragnathidae 


Metellina segmentata 


0 


1 





Prey habitat associations 

The majority of the prey MOTUs identified to the species 
level were habitat generalists. Generalist prey species tended 
to be consumed by both bat species and formed the majority 
of their diet (P austriacus 65%, P. auritus 68%). Of the 25 
prey MOTUs with specific habitat associations, the majority 
of species that the literature considered as found exclusively 
in woodlands (67%) were consumed by P. auritus, while the 
majority of species associated with more open grassland habi- 
tats (75%) were consumed by P. austriacus. Patterns of prey 
selection by habitat were not random (Chi square: x 2 = 13.1, 
df= 1,P<0.01). 

Concordance between morphological and 
molecular analysis 

A morphological analysis of 170 P. austriacus fecal samples 
resulted in the identification of six prey orders, the most com- 
mon of which was Lepidoptera (62%), followed by Diptera 
(32%). The morphological analysis identified two prey or- 
ders that were not recorded in the molecular diet analysis, 
Coleoptera and Hymenoptera, but instead did not identify 
any Trichoptera or Isopoda (Fig. 2). Nevertheless, differences 



in the proportion of prey orders were not statistically signifi- 
cant (Paired t -tests, Arcsin transformation: P = 0.98) . Similar 
to the results of the molecular diet analysis, the proportion 
of Lepidoptera in the diet differed between the study sea- 
sons, being significantly lower in spring (Kruskal-Wallis test 
with Arcsin transformation: H = 47.7, df = 2, P < 0.001), 
while the proportion of Diptera was highest in the spring 
(Fig. 4). 

A comparison with a previous morphological study of 
the diet of P auritus (data adapted from Hollyfield 1993, 
N = 240) showed similar proportions of Lepidoptera (68%) 
and Diptera (20%) but higher proportions of Coleoptera 
(6%) than those identified in our molecular study (Fig. 2). 
The morphological analysis revealed four additional prey 
orders absent from our study, Neuroptera, Trichoptera, 
Hymenoptera, and Plecoptera, though all contributed 
together to less than 1.6% of the diet. As in the diet of P. 
austriacus, there were no significant differences in the results 
of the molecular and microscopic diet analysis (P = 0.59). 
Niche overlap between the diet of P austriacus and P auritus 
at the ordinal level, based on the result of the morpholog- 
ical diet studies, was very high and significant (Ojk = 0.95, 
P = 0.01). 
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Lepidoptera 

Diptera 

Neuroptera 

Hemiptera 

Trichoptera 

Isopoda 

Coleoptera 

Arachnids 

Hymenoptera 



(a) P. austriacus - molecular study 



(b) P. auritus - molecular study 




(c) P. austriacus - morphological study (d) P. auritus - morphological study 




Figure 2. Pie charts showing the diet composition at the ordinal level of the two bat species in southern England based on molecular ([A] Plecotus 
austriacus [N = 28]; [B] P. auritus [N = 24]) and morphological diet analyses ([C] P. austriacus from the Devon colony [W = 1 70]; [D] P. auritus adapted 
from Hollyfield 1993 [W = 240]). The proportion of prey orders in the diet is presented as percent frequency of occurrence (%FO). 



Discussion 

Use of high-throughput DNA sequencing 
to study interspecific interactions 

High-throughput DNA sequencing technology has the po- 
tential to increase the scope of dietary studies by revealing 



patterns of resource partitioning among competing species 
with a level of detail not previously possible (Valentini et al. 
2009a). This study is the first application of molecular tech- 
niques to study interspecific competition between sympatric 
cryptic species that are expected to show high dietary over- 
lap due to strong morphological similarities (Findley and 
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Figure 3. The proportion of Lepidoptera families identified in the diet 
of P. austriacus (black bars, N = 28) and P. auritus (gray bars, N = 24) in 
southern England. 



Black 1983). Unlike previous DNA sequencing techniques, 
which relied on either cloning (e.g., Deagle et al. 2005; Zeale 
et al. 2011) or the selection of prey fragments (Clare et al. 
2009, 201 lb), high-throughput DNA sequencing exploits en- 
tire samples and does not rely on subsampling, and therefore 
even very rare dietary items can potentially be detected (Dea- 
gle et al. 2009). Moreover, the use of uniquely tagged primers 
(Binladen et al. 2007) makes HTS more economical by en- 
abling the pooling of many samples (up to 30 fecal samples 
per sequencing lane in our study) together while retaining 
the ability to ascribe prey sequences to an individual bat con- 
sumer, thus allowing the identification of seasonal and sexual 
dietary and niche overlap patterns. Yet, due to differences in 
primer binding, prey digestibility, and amount of DNA in 
prey tissue, neither method of molecular diet analysis to date 
can provide an accurate quantitative measure of diet compo- 
sition within any one sample (Sipos et al. 2007; Deagle et al. 
2010). 

The number of potential prey items (MOTUs) identified in 
the diet of the two bat species in our study (P. austriacus [90] , 
P. auritus [68]) did not exceed those identified in previous 
molecular studies of the diet of insectivorous bats (Clare et al. 
2009 — 127 prey species in the diet of Lasiurus borealis; Zeale 
201 1 — 89 prey species in Barbastella barbastellus). However, 
the number of MOTUs per fecal sample in our study was 
much higher (we identified a maximum of 17 prey MOTUs 
per dropping in P. austriacus and 16 in P. auritus, versus 
a maximum of seven and nine prey species in Clare et al. 
(2009) and Zeale (2011), respectively), suggesting that HTS 
has greater potential to reveal rare dietary components and 
species with less visible remains in the feces. 
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Prey orders 

Figure 4. Seasonal variation in the proportion of prey orders identified 
in the diet of P. austriacus in southern England based on (A) the molec- 
ular study (Spring: N = 8, Summer: N = 14, Autumn: N = 6) and (B) 
the morphological study (Spring: N = 60, Summer: W = 80, Autumn: 
N = 30). Diet composition is presented as percent frequency of occur- 
rence (%FO). 



The application of molecular techniques did not result 
in the identification of more lepidopteran prey species than 
previous studies of culled prey remains. Bauerova (1982) 
recorded 140 species of Lepidoptera in the diet of P. aus- 
triacus in central Europe, while Robinson (1990) found 34 
species of Lepidoptera in the diet of P. auritus in England. 
However, because DNA sequencing techniques are less likely 
to be biased toward large hard-bodied prey (Whitaker et al. 
2009), we were able to successfully record a variety of smaller 
and softer bodied prey items of orders overlooked in previous 
studies. 
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Figure 5. Extent of between (interspecific) and within species niche 
overlap among P. austriacus and P. auritus in southern England based 
on Pianka's measure of niche overlap, including the significance of niche 
overlap relative to random simulations (**P < 0.001; *P < 0.01). Within 
species niche overlap was tested by comparing the diet of Devon and Isle 
of Wight colonies separately for each bat species (P. austriacus: Devon 
N = 1 9, Isle of Wight N = 9; P. auritus: N = 1 0, 7). 

Strong overall agreement between the diet composition 
and seasonal patterns of diet composition of the two bat 
species identified in our HTS study and our own and pre- 
vious (Hollyfield [1993] for P. auritus) morphological diet 
studies supports the validity of the molecular diet composi- 
tion estimations at the ordinal level. Moreover, the fact that 
similar results were obtained from morphological studies car- 
ried out on > 170 fecal samples suggests that our sample size 
was sufficient for representing differences in the diets of the 
two bat species. 

The short sequence length (157 bp) and lack of complete 
overlap between the coverage of the primers used in our study 
and available COI reference sequences in BOLD (only ~ 130- 
bp overlap) meant that not all MOTUs could be identified 
to the species level, and we were not able to use phyloge- 
netic trees to confirm identification. In addition, many of the 
sequences showed 100% match to more than one species 
distributed globally, some of which belonged to different 
genera, thus forcing us to create a three-tiered confidence 
system and to identify some species based on whether they 
are found in the United Kingdom. Nevertheless, because a 
high proportion of prey DNA fragments recovered from fe- 
ces are degraded and relatively short (Deagle et al. 2006), 
short amplicons may provide a more accurate estimate of 
diet composition, including also prey of higher digestibility, 
and are more likely to identify rare dietary components. As 
such, short amplicons can overcome problems encountered 
by Clare etal. (2011b) of low amplification success and high 
contamination by nonprey DNA. The source material may 
thus dictate the choice of primer length, a trade-off between 
length of amplicon for identification and the impact of DNA 
degradation. 
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As suggested by Zeale et al. (201 1), the primers successfully 
amplified DNA from a wide range of arthropod orders, in- 
cluding a crustacean order (Isopoda) never before recorded 
in the diets of the two bat species. Zeale et al. (2011) fed an 
individual P. auritus an experimentally manipulated diet of 
Diptera, Coleoptera, and Lepidoptera, and found close agree- 
ment between items recorded in the diet by molecular analysis 
and the food items fed to the bat 6-24 h previously, suggest- 
ing that the primers do not appear to suffer from significant 
amplification biases. Microscopic analysis confirmed the diet 
composition results concluded from genetic data. Although 
previous morphological studies found higher proportions of 
Coleoptera in the diet of both bat species, our microscopic 
analysis of genetically analyzed fecal samples confirmed the 
low proportion of Coleoptera in our samples and the absence 
of other prey orders. Hence, the primers appear to provide 
good coverage of the main prey orders commonly found in 
the diet of insectivorous bats. 

Molecular approaches were able to overcome many of the 
biases associated with traditional bat diet analysis techniques 
(reviewed in Whitaker et al. 2009), while providing an accu- 
rate estimate of diet composition and more fine-scale extent 
of niche overlap (see Table 2 for comparison of the perfor- 
mance of the various diet analysis methods). Yet, molecular 
diet analysis is still constrained by an incomplete reference 
sequences database, though this problem was at least partly 
overcome in our study through the use of MOTUs as a surro- 
gate for prey species identification. High correspondence be- 
tween MOTUs and lepidopteran species, for which the BOLD 
reference dataset is relatively complete, suggests that MOTUs 
are an effective method of diet assessment. However, our re- 
sults suggest that they may overestimate prey species richness 
by 12% due to the proportion of repeated species identifica- 
tion in our study. This overestimation can be attributed to 
intraspecific polymorphism in the amplified region (Valen- 
tini et al. 2009b), though it can also be an artifact of the 
incomplete reference dataset or taxonomic ambiguity among 
some prey species. 

Trophic ecology relative to 
resource availability 

Insectivorous bats are sometimes perceived as opportunistic 
foragers (Kunz 1974), consuming prey based on their avail- 
ability both in the environment (Swift et al. 1985) and to the 
bat itself, given the constraints of its sensory and flight abili- 
ties (Fenton 1990; Siemers and Schnitzler 2004). Clare et al. 
(2011b) shows that the same predator can be regarded as a 
generalist, based on the number of prey MOTUs identified in 
the diet, and as a specialist when the diet is considered phylo- 
genetically, if all the prey species identified in the diet belong 
to the same order. In our study, diet selectivity at the ordinal, 
and even family level, appears to be high. Lepidoptera, and 
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Table 2. Comparison of the strength and weaknesses of available diet analysis techniques. 







Molecular techniques 


Traditional methods 




High-throughput 


Cloning (Zeale et al. 


Fragments (Clare et al. 


Morphological — 






sequencing (HTS) 


2011) 


2011b) 


feces/stomach 


Culled prey remains 


Diet resolution 


High: species level 


High: species level 


High: species level 


Low: family order 


... 

High: species order 




constrained by 


lui ibLi dii it;u uy 


lui ibLi dii ifciu uy 




\~\\ it Hononrlc r\r\ 
UUL UtptllUb UN 




reference sequences 


i cicifcri ilc btqufcfi iLtb 


rafo rp n com icin/-pc 
icIfcrlcllLc bfc!qUfc;i lLfclb 




taxa 




database (can be 


database (can be 


database {can be 








mediated by MOTU) 


mediated uy iviu I u) 


m/-i,-l i -,+nrl l~ii ( l\ /\(~\T\ \\ 

mediated uy iviu I u) 






Diet coverage 


Better representation 


Limited by the 


Selected fragments. 


Greater bias toward 


□ ias toward large prey 




of the DNA 


selected number of 


Bias toward prey of 


hard-bodied prey of 


that requires culling 




extracted, but bias 


i"li*"\nac Riac tr"UA/arrl 
LIUIIfcfb. Dldb LUVVdlU 


Ii"\\a/ rl inpct ihi 1 i+v nr 

iuw uiy trbLiuiii Ly ui 


iuw uiy tbiiuiii Ly 






from primer binding 


f mo f lf~i a + i nn Pi M A 
Wcc IIUdLlliy UINM, 


f m an + cal orti r~ir\ 
lldyillfcrllL btrlcLLIUri 








biases, and relative 


DNA abundance, 


method 








DNA abundance of 


dllU pillllcl Uldb 










prey 










Diet quantification 


Number of species/ 


Number of species/ 


Number of species/ 


Number of prey taxa 


Number of prey taxa 




MOTUs, no within 


iviu i us, no wi li nil 


iviu i us, no w i Li 1 1 r 1 


dllU peiLfcfllL VUlUllltf 


dMU Ufcrl Ibl Ly 




samples 


c^m nipt: 
bdl 1 1 [Jlcb 


ca m nipt; 
bdl 1 1 [Jltrb 








quantification 


quantification 


quantification 






Rare dietary 


High potential to 


Selection of clones for 


Potential of 


High potential to 


Rare prey identified 


components 


identify rare 


sequencing reduces 


identification 


identify 


only if require 




prey — high 


the chance of 


depends on number 


hard-bodied rare 


culling 




proportion of DNA 


lueruiiyiny iditr prey 


Ol 1 1 dyi lltf II Lb 


piyy, UUL IOW 






variance sequenced 




sequenced from 
each dropping 
(costs) 


potential to recover 
idie boi l picy 




Accuracy 


Does not require 


U'Ucj iiul rtruuiitr 


rinoc nnt romnro 
l-'Utrj MUL Itruuiltr 


Ar~r"i irar\/ roni nrpc 
MLLUIdLy ItrLjUllcb 


r\<_LUIa<_y ItrLIUlrtrj 




taxonomists to 


tavnnn m i c1"c ir\ 
LdXUI IUI 1 1 lb Lb LU 


t^vnnnm i etc ir\ 
LdAUl IUI 1 1 lb Lb LU 


cinn if ii-an+ 
biy 1 III ILdl 1 L 


fcrl 1 LUI 1 lUIUy ILdl 




obtain accurate 


obtain accurate 


OULdlll dLLUrdLfcf 


tillLOIllOIUyiLdl 


traminn r\r\ + K ja r^a r+ 
Lldllllliy Oil Lllfcf pdl L 




results, lower 


results, lower 


results, lower 


training on the part 


of the identifier 




analyzer bias 


analyzer bias 


analyzer bias 


of the identifier 




Costs 


Less limited by the 


Expensive cloning 


Depends on number 


Very low — minimal 


Very low— minimal 




amount of prey 


process, but less 


of fragments 


consumables 


consumables 




sequences or 


specialized facilities 


sequenced — each 








droppings per 




fragment requires 








sequencing run. 




separate 








LUWcl LUbL UfcM 




sequencing. Less 








sequence but high 




specialized facilities. 








initial costs and 












needs specialized 












facilities 










Applicability across 


High relevance across 


High relevance across 


Mainly relevant for 


Mainly relevant for 


Only relevant for 


feeding groups 


taxa, including 


taxa, including 


predators and seed 


predators 


predators that use 




herbivores 


herbivores 


disperser 




feeding perches 



in particular eared-moths of the family Noctuidae, were the 
main prey consumed by both bats. In contrast, when the diet 
is considered at the MOTU or species level, P. austriacus and 
P. auritus appear to have a generalist diet, including a variety 
of prey items, most of which were often only consumed by 
a single individual (although if more samples could be ob- 
tained more individuals may have consumed these items). 
Diet generality is reinforced by the fact that all prey MOTUs 
identified to the species level were common and widespread 



in southern England, and were consumed by the bats when 
at the peak of their adult-form activity (Heath and Emmet 
1983; Waring etal. 2003). 

The larger size and potentially high energetic value of Lep- 
idoptera makes them an important component of the diet 
of many insectivorous bat species (e.g., Vaughan 1997; Bog- 
danowicz et al. 1999). Based on a small sample size, we found 
a trend in the diet of P. austriacus and P. auritus to shift 
from a diet dominated by Lepidoptera in summer toward the 
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consumption of high proportions of Diptera during spring 
and autumn, when Lepidoptera availability is reduced (Jones 
1990). This seasonal trend of dietary shift highlights the pref- 
erential consumption of Lepidoptera and the tendency of P. 
austriacus and P. auritus to concentrate on high-quality prey 
at times of high energetic demands (Kurta et al. 1989; McLean 
and Speakman 1999, for P. auritus). 

Many moths have tympanal organs that enable them to 
detect the echolocation calls of approaching bats and initi- 
ate behavioral responses to evade predation. They are par- 
ticularly sensitive to relatively high-intensity echolocation 
calls between 20 and 60 kHz (Rydell et al. 1995; Miller 
and Surlykke 2001). Consequently, the proportion of Lep- 
idoptera is higher in the diet of bats that either produce 
echolocation calls above or below these frequencies (Bog- 
danowicz et al. 1999), or produce low-intensity echolocation 
calls when approaching prey (Goerlitz et al. 2010), or listen 
to prey-generated sounds to detect prey (Faure and Barclay 
1992). High proportions of Lepidoptera in the diets of P. aus- 
triacus and P. auritus have been attributed to the latter two 
strategies (Rydell et al. 1995) due to the bats' low-intensity 
echolocation calls, long ears, and the ability of P. auritus to 
glean prey from surfaces without producing echolocation 
calls (Coles et al. 1989; Anderson and Racey 1991). How- 
ever, P. austriacus appears to hunt mainly by aerial hawk- 
ing rather than gleaning (Bauerova 1982) and the use of 
passive listening by this species is yet to be confirmed (Swift 
1998). 

Mechanisms of coexistence in cryptic 
bat species 

Although dietary partitioning is regarded as a principal mech- 
anism for resource partitioning in animal communities (e.g., 
Schoener 1974; Saunders and Barclay 1992, for bats), it does 
not appear to play an important role in facilitating coexistence 
among sympatric P. austriacus and P. auritus. High similarity 
in both morphology (Ashrafi et al. 2010) and echolocation 
call characteristics (Russo and Jones 2002) suggests that the 
two species may not differ sufficiently in their aerodynamic 
or sensory abilities to have access to different prey (Siemers 
and Schnitzler 2004). Significant dietary overlap even when 
prey resources were limiting led Schoeman and Jacobs (201 1) 
to conclude that competition does not structure the trophic 
niches of bat species, although they admit that greater 
differences in diet may be apparent when diets are quan- 
tified at finer resolution than at the ordinal level. Indeed finer 
resolution of dietary composition revealed the importance of 
the partitioning of minor prey items in reducing interspecific 
competition among coral reef fish, despite high overlap in the 
consumption of major prey items (Nagelkerken et al. 2009). 
Similarly, high dietary overlap in our study was primarily the 
result of the sharing of common prey items between the two 
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cryptic bat species, while rare dietary items were generally 
consumed by a single species (Table 1). 

We found that interspecific prey resource overlap was 
higher than intraspecific overlap, suggesting that individual 
P. auritus bats, in particular, are more likely to consume the 
same prey items as individuals from their sibling species, P. 
austriacus, than members of their own species. Greater lev- 
els of intraspecific dietary partitioning may be the result of 
higher intraspecific than interspecific prey resource compe- 
tition because of differences in foraging habitat use. While 
P. austriacus preferentially forages in more open habitats, 
such as unimproved grasslands, or at the edge of riparian 
areas and woodlands (Razgour et al. 2011), P. auritus forages 
primarily in deciduous woodland (Entwistle et al. 1996). In 
this case, apparent resource partitioning may be less of an 
effect of competition for limited food resources than a pat- 
tern that results from habitat partitioning. The majority of 
grassland-only prey species, such as the moths Agriphila tris- 
tella, Chrysoteuchia culmella, and the crane fly Tholera deci- 
malis, were consumed by P. austriacus, while woodland-only 
moth species, including A. epomidiona.ndAcronictaalni,weie 
consumed by P. auritus. Generalist prey species, such as N. 
pronuba, the most common prey item in the diet of both bat 
species, tended to be consumed equally by both species, hence 
resulting in the observed high dietary overlap. Other cryptic 
Myotts and Pipistrellus bat species show similar patterns of 
foraging habitat segregation mirrored by prey habitat asso- 
ciations (Arlettaz 1999; Davidson-Watts et al. 2006), though 
in both cases dietary partitioning was significant despite the 
coarse study resolution (Arlettaz et al. 1997; Barlow 1997). 

The diet of the two cryptic bat species tended to only over- 
lapped in summer, when the availability of Lepidoptera, the 
preferred prey resource, peaks in southern England (Jones 
1990), suggesting that dietary resource partitioning is great- 
est in autumn, when prey resources are more limiting. In 
addition, interspecific competition may be greater in au- 
tumn because P. austriacus tends to forage in woodlands to a 
greater extent when ambient temperatures are low (Razgour 
et al. 2011). Therefore, it appears that when prey resources 
are more limiting and habitat use overlap is higher, the role of 
dietary partitioning in facilitating coexistence becomes more 
important. However, a larger dataset is needed to confirm the 
trend of seasonal differences in dietary niche overlap identi- 
fied in our study. 

Our results suggest that instead of prey selection driv- 
ing differential habitat use (as is the case in Saunders and 
Barclay 1992), sympatric P. austriacus and P. auritus sim- 
ply consume the most energetically beneficial prey available 
to them in their respective foraging habitats. Hence, spatial 
partitioning, rather than dietary partitioning, is the princi- 
pal mechanism of coexistence among these highly cryptic 
bat species. It is still unclear what morphological differences, 
if any, drive the differences in foraging and flight behavior, 
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because the two bats show no significant differences in wing 
loading and aspect ratio (Sevcik 2003), the main wing pa- 
rameters controlling manoeuvrability and ability to sustain 
prolonged flight in bats (Norberg and Rayner 1987). Indeed, 
it could be that morphological divergence is not necessary to 
drive differences in habitat use by these species. 

Conclusions 

HTS provided species-level identification of dietary compo- 
nents, and therefore allowed a more comprehensive compari- 
son of the diet of two sympatric cryptic bat species. We found 
little evidence of trophic resource partitioning. Instead, pat- 
terns of prey consumption corresponded to predator habitat 
use, highlighting the importance of habitat partitioning as 
the principal mechanism of coexistence. The use of MOTUs 
as surrogates for prey species enabled us to overcome the 
constraints of an incomplete reference dataset, and corre- 
sponded well with prey species identification. Using generic 
arthropod primers, we were able to PCR amplify and se- 
quence short fragments of degraded DNA from bat feces 
with little contamination from nonprey DNA, but our prey 
species identification confidence was compromised. There- 
fore, efforts should be put into testing and developing an 
optimal amplicon length for molecular diet studies, which 
will provide maximum identification ability with minimum 
loss of degraded prey DNA. Our study provides significant 
support for the use of molecular diet analysis, and in par- 
ticular HTS, to decipher potential mechanisms of resource 
partitioning. 
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